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Abstract 


We consider the most general set of SU{2) x U{1) invariant CP-violating operators of di¬ 
mension six, which contribute to VVh interactions {V = W,Z,^). Our aim is to constrain 
any CP-violating new physics above the electroweak scale via the effective couplings that arise 
when such physics is integrated out. For this purpose, we use, in turn, electroweak precision 
data, global fits of Higgs data at the Large Hadron Collider and the electric dipole moments 
of the neutron and the electron. We thus impose constraints mainly on two-parameter and 
three-parameter spaces. We find that the constraints from the electroweak precision data are 
the weakest. Among the existing Higgs search channels, considerable constraints come from 
the diphoton signal strength. We note that potential contribution to h ^ jZ may in principle 
be a useful constraining factor, but it can be utilized only in the high energy run. The contri¬ 
butions to electric dipole moments mostly lead to the strongest constraints, though somewhat 
fine-tuned combinations of more than one parameter with large magnitudes are allowed. We 
also discuss constraints on gauge boson trilinear couplings which depend on the parameters of 
the CP-violating operators . 
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1 Introduction 


Although the discovery of “a Higgs-like boson” at the Large Hadron Collider (LHC) has been a 
refreshing development 00 . there is no clear signal yet for physics beyond the standard model 
(SM). It is therefore natural that physicists are trying to wring the last drop out of the Higgs 
sector itself, in attempts to read fingerprints of new physics. 

One approach is to examine all available data in terms of specific new models, such as supersym¬ 
metry or just additional Higgs doublets. In the other approach, one can take a model-independent 
stance, parametrize possible modifications of the interaction terms of the Higgs with pairs of SM 
particles, and examine them in the light of the available data. Such modifications can again be of 
two types. In the first category, they are just multiplicative modifications of the coupling strengths, 
the Lorentz structures remaining the same as in the SM. Constraints on such modifications have 
already been derived from the available Higgs data [3]-[^. In the second class, one considers ad¬ 
ditional operators with new Lorentz structures satisfying all symmetries of the SM Gauge 

invariance of such operators in their original forms may be expected, since they are obtained by 
integrating out new physics that is just above the reach of the present round of experiments. Sets 
of such higher-dimensional operators contributing to the effective coupling of the Higgs to, say a 
pair of electroweak vector bosons have been studied extensively. Here it makes sense to include 
only SU{2) X U{1) invariant operators in one’s list to start with, because the yet unknown new 
physics lies at least a little above the electroweak symmetry breaking scale. A host of such gauge 
invariant higher-dimensional operators have been, and are being, analyzed with considerable rigor, 
and now there exist limits on them, using data ranging from electroweak precision measurements 
to global fits of LHC results [T6] - [4^ . 

Most of such analyses include higher-dimensional operators that conserve charge conjugation 
(C), parity (P) and time-reversal (T). However, there is evidence of C, P and CP-violation in weak 
interactions (^, and there are speculations about other sources of CP-violation as well, especially 
with a view to explaining the baryon asymmetry in our universe |^[^ . The possibility of CP- 
nonconservation cannot therefore be ruled out in the new physics currently sought after. Thus one 
may in principle also obtain higher-dimensional interaction terms involving the Higgs and a pair of 
gauge bosons. The constraints on such terms, and identification of regions in the parameter space 
where they can be phenomenologically significant, form the subject-matter of the present paper. 

The CP-violating effective couplings, interestingly, are not constrained by the oblique elec¬ 
troweak parameters at one-loop level up to O(^), where A is the cutoff scale of the effective 
theory. The leading contributions to self-energy corrections to electroweak gauge bosons at one- 
loop level occur at Therefore, electroweak precision (EWP) data are not expected to 

provide severe constraints on CP-odd parameters. In addition to this, Higgs-mediated event rates 
in various channels receive contributions from these couplings at Thus they can also be 

constrained from global fits of the LHC data. The strongest limits on them, however, arise from 
the contributions to the electric dipole moments (EDMs) of the neutron and the electron, both of 
which are severely restricted from experiments. As we shall see in the following sections, a single 
CP-violating operator taken at a time may in certain cases be limited to a very small strength 
from the above constraints, while two or three such operators considered together can have rel¬ 
atively larger, but highly correlated coefficients. Some of these operators can have interesting 
phenomenological implications, especially in the context of the LHC. 

A study in similar lines can be found in Refs. 


42,45-47,49 . However, we have performed the 


most comprehensive analysis, taking all the five possible dimension -6 CP-violating VVh operators 
{V = VE, Z, 7 ), which are not yet discussed thoroughly in the literature. We provide the constraints 
obtained from the oblique electroweak parameters. The constraints coming from global fits of LHC 
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data and electric dipole moments, for two and three operators taken at a time, have been compared. 
In addition to these we also provide the constraints on trilinear gauge boson couplings coming from 
LEP data on gauge boson pair production. 

Our paper is organized as follows. In section 2, we provide details on the CP-violating gauge 
Higgs operators and derive Feynman rules for three point vertices of our interest. In section 
3 we present the constraints on CP-violating parameters coming from the precision electroweak 
measurements. Following this, we perform a global analysis of these parameters using LHC data 
on Higgs in section 4 and we discuss EDM constraints on CP-violating parameters in section 5. 
Section 6 is devoted to discussion of results. Finally we summarize and conclude in section 7. 


2 Morphology of CP-violating Gauge-Higgs Operators 

In the effective Lagrangian approach that has been followed here, one can write a Lagrangian (>Cefr) 
comprising only of the SM fields, where the effects of the new physics that appear above the cutoff 
scale A are encapsulated in higher dimensional gauge invariant operators. In general, 

= E <i) 


where > 4 is the mass dimension of the operator and the dimensionless free parameter fi fixes 
the strength of the corresponding operator. The operators constructed out of the Higgs doublet 
and the SU{2) x U{1) gauge fields are of even dimensions, and at the leading order they have 
mass dimension di = Q. The dimension six CP-even gauge invariant operators constructed out of 
the Higgs doublet ($) and the electroweak gauge fields that modify the gauge-Higgs 

couplings are given as follows: 

Ow = Os = 

Obb = Oww = 

Obw = (2) 


In the above, we have defined and 9 d' electroweak 

coupling parameters corresponding to SU{2) and U{1) gauge groups respectively, and r^{a = 
1,2,3) are the three Pauli matrices. We define the gauge covariant derivative as = d^ — 
— where Y is the hypercharge quantum number. With this choice of the definition 

of gauge covariant derivative, field strength tensor is given by, VF^^ = d/jW^ — dyW^ + 




The constraints on CP-even parameters and their collider implications have been 


studied extensively in the literature 16,21,32,36,43,44,48 . In this work we are interested in 


corresponding CP-violating dimension six gauge-Higgs operators. These are. 


Ow = Os = 

Obb = Oww = 

Obw = (3) 
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Coupling 

Effective coupling strength 

CwWh 

{ — fw — “^fww) 

Czzh 


O'j'jh 

— + Jbb — Ibw) 

C^zh 

iw/^ { — fw + Jb) + ^{SwfBB — Cwfww) + ‘^C2wfBW 

CwW'y 

Swl‘^{fw + fs + ‘^Ibw) 

Cwwz 

-Swtw l‘^{fw + /b + ‘^fBw) 


Table 1: CP-odd VVh and WWV coupling factors and their effective strengths. 




being the four-dimensional fully 


where, and 

antisymmetric tensor with = 1. 

In principle, the CP-even operators [Eq.Q] could have been assumed to exist simultaneously 
with the CP-odd ones considered here. However, such an approach generates far too large a set 
of free parameters, where the signature of the CP-violating effective couplings would be drowned. 
Moreover, the CP-even operators are independent of the CP-odd ones (and vice versa); therefore, 
setting them to zero is a viable phenomenological approach. We therefore postulate that the new 
physics above scale A is such that only CP-violating dimension six effective operators are apprecia¬ 
ble^ and the corresponding CP-conserving ones are much smaller. Such a “simplified approach,” 
we reiterate, is unavoidable for unveiling CP-violating high scale physics, as has been recognized 
in the literature (49,53-^. Studies focusing exclusively on the generation of CP-violating terms 


in specific new physics frameworks can also be found, an example being those in the context of 
extra space-time dimensions . 

Since we focus on the extension of the SM through the inclusion of the CP-odd operators only, 
the full BSM Lagrangian looks like, 

^BSM = >CsM + Ow + Oww + Ob + Obb + Obw, (4) 

where CsM is the standard model Lagrangian. After the electroweak symmetry breaking, these 
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CP-odd operators contribute to following three-point vertices of our interest^ 

gmw, 


^WWh 

^ZZh 


L 


77/1 


L. 


-iZh 


A2 

gmw 

A 2 


ifw + 2 fwwV’''^^kiak 2 ^W+{h)W-{k 2 )h{k), 
C^fw + 'S^/b , ‘^{c'^fww + S'^fsB) 


( 5 ) 


-b 






-W 


e^^^^kiak2^Z^(ki)Z,(k2)h(k), 

')^w(fww + Jbb — Ibw) 

(“/vK + Ib) + ‘^{s^rfsB - c\rfww) + ‘^C 2 wfBW 


lww-1 = ( ^ 2 A^ + /b + ‘^Jbw) 

(fci) fp- ik2)Aa{k), 


( 6 ) 


(7) 


( 8 ) 


(9) 


k^wwz 


= ~( ^2A^ ^iswtw){fw + /b + “^Ibw) 
e^^’^^^kpW+ (fci)fp- (fe2)^a (fc). 


( 10 ) 


In the above equations sw = sin^vv? = cos^vv? tw — tan^vv? 2 ind C 2 w — cos29w^ where 9w 
is the Weinberg angle. Here ks are the four-momenta of the fields that enter the vertex. We 
have taken all momenta to be inflowing toward the three-point vertex in establishing the Feynman 
rules. From the list of CP-odd interaction vertices shown above, one can observe a general tensor 
structure of the form e^^^^kiak 2 ^ in VVh vertices and a general tensor structure of the form 
in trilinear gauge boson couplings (WWV). Because of this the CP-odd couplings are 
linear combinations of the parameters /^. Note that we have not included the CP-odd operator 
involving gluon-Higgs coupling, 

Ogg = ( 11 ) 


A2 


This operator introduces a 9qcd term 1^,59 , and it is severely constrained by the experimental 


measurement of neutron EDM . In Tablej^ we list various couplings and their effective strengths 

ignoring the overall dimension full factor of in Cyyh and the dimensionless factor of 
in Cwwv couplings. Note that only ZZh and yZ/i couplings receive contribution from all five 
CP-odd operators. The operators which contribute to WWj also contribute to WWZ and these 
couplings are related by, Cwwz — —twCwW'y 


3 Constraints from Electroweak Precision (EWP) Data 

We note that unlike some of the CP-even (D=6) operators, the CP-odd operators do not contribute 
to the gauge boson propagator corrections at tree level, hence are not expected to receive severe 

^ The CP-odd operators considered here also contribute to four-point and five-point vertices like VVhh, VVVh 
and VVVhh. However, as we will see in following sections, all the observables used in our analysis are sensitive to 
only three-point vertices at the leading order. 
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Figure 1: One-loop self energy corrections to electroweak vector bosons (oblique corrections) in 
presence of CP-odd operators. The blobs show the effective CP-odd vertices. These corrections 
are of (9(1/A^). 


bounds from the electroweak precision data. This is due to the antisymmetry of the epsilon 
tensor which is present in all CP-odd operators. In fact, because of the same reason, all quantum 
corrections to gauge boson two-point functions up to O(^) vanishand first nonzero contributions 
due to CP-odd operators appear at (9(^). As we will see in section 4, the CP-odd couplings 
contribute to observables related to LHC Higgs data at this order. It would be interesting to 
discuss the implications of the electroweak precision measurement constraints on the parameters 
of CP-odd operators. 

It is well known that the dominant effects of new physics can be conveniently parametrized in 
terms of Peskin-Takeuchi parameters [^. These are related to the gauge boson two-point functions 
as. 


aS 

= 444(n'^^(0) - ny(0) - 

V Sw^w ^ 

(12) 

aT 

2 2 

(13) 



aU 

= ~ ^wn^^(o) — s^n(^^(o) — 2cu/su/n(^^(o)^ 

(14) 


where, IlviV 2 {p‘^) ^tnd are the part of the two-point function and its derivative with 

respect to respectively. The relevant one-loop Feynman diagrams are shown in Fig. [ij We 
have regularized ultraviolet (UV) singularities of these diagrams in dimensional regularization 
(DR). The expressions for IIy^v 2 {p^) terms of standard one-loop scalar functions are given in 
Appendix]^ We find that n^7(0) = n7^(0) = 0 which is expected due to the transverse nature of 
the photon. The renormalization of UV singularity is carried out in MS scheme which introduces 
scale dependence in these expressions. We have identified the renormalization scale with the cutoff 
scale A. Thus the gauge boson two-point functions also have ln(A) dependence apart from the 
overall 1/A^ dependence coming from CP-odd couplings. Because of this an explicit choice of A is 
necessary in deriving the EWP constraints on /^s. 

^These corrections are proportional to (p being the four-momentum of gauge boson) which is zero. 
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Figure 2: Constraints from electroweak precision data keeping two parameters nonzero at a time 
and for A = 1 TeV. 
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For A = 1 TeV, the Peskin-Takeuchi parameters due to CP-odd couplings are given by, 


5 = (-3.36 - 1.28 C^^zh + 4.64 Clzh “ 4-49 C^^hC^zh - 6.21 C^zhCzzh) x 10“^ (15) 

T = —9.74 X 10 ^ (16) 

U = (—0.960 — 4.69 C‘^zh ~ 4.64 C'zzh + 6-67 + 2.76 

— 3.59 C^^hC^zh — 4.96 CjZhCzzh) x 10“^. (17) 

We can also express them in terms of /jS using their relation with CjS given in Table The exper- 
imental limits on T and U parameters are obtained by fitting the data on various electroweak 
observables with these parameters. The limits are [^ , 


S = -0.03 ± 0.10, T = 0.01 ± 0 . 12 , U = 0.05 ± 0.10. 


(18) 


In Eqs. ( 15), (16) and ( [T7| the coefficients of various couplings are ^ 10“^ suggesting that the 
EWP constraints cannot be very strong. Therefore, we only consider the case where any two out 
of five parameters are nonzero. In Fig. we display the allowed range for CP-odd parameters 
which satisfy the above limits on T and U parameters for all ten sets of two parameters taken 
together. Here we have varied parameters freely to ensure that we obtain a bounded region. We 
note that large values of 0(1000) for /bb are always allowed. On the other hand the allowed range 
for fww never goes beyond 250. Allowed values for all other parameters can be of 0(100 — 1000 ). 
Some of these observations can be understood once we express the T and U parameters in terms 
of fiS. As we turn on other parameters, these constraints become weaker. Also, for a larger cutoff 
scale the allowed parameter space grows as one would expect. 


4 Constraints from LHC data 


The presence of CP-odd operators introduces modifications in the strength of the gauge-Higgs 
couplings, and hence changes the Higgs production and decay rates in channels involving these 
couplings. Since we are interested in CP-even observables, the SM VVh couplings which are 
CP-even, do not interfere with the CP-odd VVh couplings. Hence the lowest order (tree level) 
modifications to the decay widths (F) and production cross sections {a) are of the order To 
quantify these changes we define the following ratios for various decay and production channels. 


F^S^(/i^y) 

FSM(/^^y) 

^ h) 
(jSM(X ^ h) 


(19) 

( 20 ) 


where Y and X are used to label the final state and initial state particles in the Higgs decay and 
production channels respectively. 


4.1 Higgs decay channels 

In the SM, the 126 GeV Higgs boson predominantly decays into bb and VFVF* followed by gg^ r^r~ , cc 
and ZZ*. It also decays to 77 , 7 Z and ii~ with much suppressed rates. Out of these, h gg^ 

/i ^ 77 and h 7 Z are loop-induced decay modes and hence are sensitive to new physics. The 
decay channels which are affected by the CP-odd operators are h 77 , h 7 Z, h VFVF* and 






h ZZ*. The expressions for the ratio of the decay widths, aij in various two body Higgs decay 
channels are as follows 0 

= l + 2.84(^) 

( 22 ) 

a^z = 1 +0.856 (^) 

(23) 

= 1 + 3.35 X 10-6 (^™) 

(24) 

= l + 3.56xl0-'=(^>™) 

(25) 

azz* = 1 + 1.40 X 10-6 (^^|f^) 

(26) 

a4^ = 1 + 1.54 X 10-6 (f|f^) 

(27) 

Since the set of gauge-Higgs operators considered in our analysis do not alter the Higgs coupling 
with gluons {g) and fermions (/), we have Oiff — agg = 1. The ratios 0 ^ 2 ^ 21 / [Eq.( 25)] and 
[Eq.( 27)] correspond to the h VFVF* ^ 2i2iy and h ZZ* ^ M respectively. Here 1 stands for 
electron and muon, and v for corresponding neutrinos. The ratios aww* stud azz* which include 
both leptonic and hadronic decays of W and Z bosons are used in calculating modified total Higgs 
decay width. As mentioned earlier, the modifications to Higgs partial decay widths at leading 
order are (9(1/A^). It is in contrast to the case of CP-even dimension six gauge-Higgs operators 
where such modifications occur at 0(1/A^). Unlike WWh and ZZh couplings, the 77/1 and 'jZh 
couplings are loop-induced in the SM. In the presence of CP-odd operators these vertices receive 
contributions at tree level. This explains the relatively large coefficients in the expressions for 


[Eq.( 1^] and ar^z [Eq.( |^] as compared to the other decay width ratios. This would imply 
most stringent constraints on the parameters contributing to these decay channels. For further 
discussion on CP-odd vs. CP-even operators, we refer the reader to section 6 . 

4.2 Higgs production channels 

At the LHC, the dominant mode to produce Higgs boson is gluon-gluon fusion (GGF) mediated 
by a top quark loop. The other major production channels include: vector boson fusion (VBF), 
associated production with a weak boson (F/i) and associated production with a pair of top quark 
{ilK). Except GGF and tth production channels, all other channels are affected in presence of 
anomalous gauge-Higgs CP-odd vertices. Like the decay width ratios, the production cross section 
ratios also receive modifications at (9(1/A^). The ratios of the Higgs production cross sections, 
in various channels at \fS = 8(7) TeV LHC are given below. 

^ We disagree with the CP-odd part of the analytic expression for a^z in Eq. (3.17) of Ref 
expression for CP-odd term in the notations of Ref. turns out to be, 

I A\/2'K‘^a2 12 

\ + Aw)\ 

where a 2 can be identified with the factor \CryZh in our notation. 


45 


. The correct 


( 21 ) 
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Ipp^Wh 


( 28 ) 


= 1 + 5.61(5.37) X 10"^ 
jpp^wh^hiu = 1 + 5.67(5.16) X 10-4 (^|p) (29) 

jpp^Zh = 1 + 4.09(3.92) X 10-4 + 2.45(2.32) x 10-4 

+ 2.55(2.44) X 10-4 (^?) 

-P‘2 -p2 -p2 

7VBF = 1 + 7.02(5.62) x 10-® + 1.50(1.44) x 10-4 (M) + 1.84(1.80) x 10-^ 

+ 6.98(6.75) X 10-4 +4.39(4.38) x 10-*^ -1.32(1.14) x 10-^ 

+ 8.96(22.2) X 10-7^ -faJ^^ j -5.06(5.03) x 10-^ 

+ 2.63(2.59) X + 8.98(9.21) x 10"’^ 

+ 6.22(5.98) X 10-4_ 2.46(2.71) x iq-5 

- 3.51(3.60) X iq-5 +3.27(3.85) x iq- 5f^ fBBfww -^ 

- 1.43(1.45) xl0-4( ^^^{^^ ) (31) 

Ipp^tth = 7GGF = 1- (32) 

These expressions have been obtained by computing the SM and BSM cross sections at tree level 
using Madgraph under the assumption that the K-factors (due to higher order corrections) 
are same in the SM and BSM cases. For that we have implemented our effective Lagrangian in 
FeynRules and used the generated UFO model file in Madgraph. The cross sections have 
been calculated using cteqGll parton distribution functions and with default settings for 
renormalization and factorization scales. 

We would like to point out that there is an additional diagram which contributes to pp ^ Zh 
due to tree level CP-odd 7 Z /1 coupling. Similarly, in VBF channel additional diagrams appear 
due to both 77/1 and 7 Z /1 couplings. Because of a different parametrization, this information is 
not explicit in the expression for VBF. We find this parametrization more convenient in terms of 
evaluating the coefficients in Eq.(|^. Also, the VBF coefficients reported above do not have any 
Vh contamination and this can be ensured in Madgraph at the process generation level. One can 
notice that the modifications induced by the CP-odd operators are relatively weak because the SM 
cross sections are already tree level effects in the modified production channels. 

One important fact in relation to these production cross section ratios ( 7 s) is that the numerical 
coefficients present in these expressions are very much cut dependent. This is associated with the 
fact that the anomalous couplings induced by the gauge-Higgs operators have a different Lorentz 
structure (therefore, different kinematic dependence) than their SM counterparts. Hence in general 
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Production channel 

8 TeV cross section (pb) 

7 TeV cross section (pb) 

GGF 

18.97 

14.89 

VBF 

1.568 

1.211 

Wh 

0.686 

0.563 

Zh 

0.405 

0.327 

tth 

0.1262 

0.0843 

bbh 

0.198 

0.152 


Table 2: Higgs production cross section in the SM 


64 


the SM and BSM cut efficiencies are not the same for a given process. The differences between 
the two become more pronounced for higher values of the CP-odd couplings. But for reasonably 
low values of the same one can still work under the approximation that the two cut efficiencies 
are the same. In this work we have taken this approximation into consideration, and taken only 
default cuts in Madgraph to simulate any production or decay channel. Since we have taken data 
only from individual production channels and not from combined channel data (e.g. Vh combined 
channel) in our chi-square analysis, this approximation finds a stronger footing. 

4.3 Global analysis 

The quantitative measure of the difference between the Higgs data from the LHC, and its corre¬ 
sponding SM predictions is given by what we call the signal strength, defined as, 

^ h)BR^^^{h ^ F) 

^ ^ h)BR^^{h^Y) ’ 

where, BR(/i Y) = is the branching ratio for Higgs decaying into Y final state, and 

^totai is the total Higgs decay width. For 126 GeV Higgs boson F^^^^ = 4.2 MeV. The total Higgs 
decay width in the BSM construct F^^^ can be expressed in terms of the SM total Higgs decay 
width rSM / by, 

■pBSM Q -pSM /q/IA 

^ total — ^total^ total' 
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Stotal is given in terms of the various branching fractions of the Higgs in the SM as, 


Stotal 


opSM I opSM 


+ BR™ + + a^zBRl 


+ 0!ww*^^n^* + 0!ZZ*^^Wz* 

which becomes on solving, 

Stotal ^ 0.736 + 0.0023(a^-^ + O.OOlGcr^^ + 0.23c^ww* + 0.029(yzz* • 


"^77 

>SM 


SM 


( 35 ) 


(36) 


The SM branching fractions for 126 GeV Higgs are taken from 64 . The signal strength in Eq.([ 
can be rewritten in a compact form using the decay and the productions cross section ratios defined 
above. 


= 7x 


ay 

Stotal 


(37) 


To perform the global fit of our CP-odd parameters, we use the standard definition of the 
chi-square function. 




- 1 . 


V' L 


^x,r 


(38) 


where fif^ is the theoretical signal strength expected in presence of CP-odd operators, and 
is the experimental signal strength reported by the LHC experiments. is the experimental 

uncertainty in fjtexp • The experimental data reported generally has unsymmetrical uncertainties 
'^XY '^XY' that we use symmetrizes these uncertainties through the following 

definition. 


= 



(39) 


Since the LHC data that we use includes data from both 7 and 8 TeV LHC runs, the theoretical 
signal strength in Eq.( [SS]) is obtained after combining the signal strengths calculated for 7 and 8 
TeV LHC. For that we have used following formula [43] , 


l^th 








a 




(40) 


where C 7 and Cs are the luminosities at 7 and 8 TeV, respectively, and and are the SM 
cross sections at those energies. These cross sections are listed in Table 

In our analysis we have taken total 15 data points which are the most updated ones. We have 
listed them in Table Note that due to large uncertainty we do not include Higgs data in h ^ 7 Z 
decay channel from CMS and ATLAS. The global analysis with five CP-odd parameters 

results into Xmin ~ 6.78. However, in this case the best fit point is very unstable with respect to 
the step size that we choose to scan the parameter space. Also, we find that the Xmin insensitive 
to the parameters fs and fw- These are the parameters which do not enter the 77/1 vertex and 
their coefficients are very small compared to those of /bb, fww and fsw^ which do enter the 77/1 
vertex [Eq.([^]. Since the LHC observables have an overall cutoff scale dependence, the ratio 
/i/A^ can be taken as the effective parameter to be constrained. In other words, the constraints 
from global analysis can be easily predicted for any value of A of interest. 

We organize the constraints on CP-odd parameters from global fit of LHC data in the following 
two parts. 
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Production channel 

Decay channel 

Signal strength 

Energy in TeV 
(Luminosity in fb~^) 

GGF (ATLAS) 

/i ^ 77 

1.32 ±0.38 6 

7 

7(4.5) ± 8(20.3) 

VBF (ATLAS) 

/i ^ 77 

0.8 ±0.7 

67 


7(4.5) ± 8(20.3) 

Wh (ATLAS) 

/i ^ 77 

1.0 ± 1.6 

67 


7(4.5) ± 8(20.3) 

Zh (ATLAS) 

/i ^ 77 

O-l^ai 67 


7(4.5) ± 8(20.3) 

tih (ATLAS) 

/i ^ 77 

1 -61 l8 67 


7(4.5) ± 8(20.3) 

GGF (CMS) 

/i ^ 77 

l-127:i 68 


7(5.1) ± 8(19.7) 

VBF (CMS) 

/i ^ 77 

1.587:6^^ 68 


7(5.1) ± 8(19.7) 

tih (CMS) 

/i ^ 77 

2.69711 68 


7(5.1) ± 8(19.7) 

GGF + tih + bbh (ATLAS) 

h ZZ* Ai 

l-7la4 69 


7(4.5) ± 8(20.3) 

GGF + tih (GMS) 

h ZZ* Ai 

r| O+0.46 

70 


7(5.1) ± 8 (19.7) 

GGF (ATLAS) 

h WW* 2l2v 

1 01+0-27 71 


7(4.5) ± 8(20.3) 

VBF (ATLAS) 

h WW* 2l2v 

1.287^3 71 


7(4.5) ± 8(20.3) 

GGF (CMS) 

h WW* 2l2v 

0.747|2 72 


7(4.9) ± 8 (19.4) 

VBF (CMS) 

h WW* 2£2u 

0 6+°-5^ 

72 


7(4.9) ± 8 (19.4) 

Wh ^ h£u (CMS) 

h WW* 2£2u 

O.bQtlil 72 


7(4.9) ± 8 (19.4) 


Table 3: LHC data used in the global analysis. 
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Figure 3: Global fits of the CP-odd parameters keeping two parameters nonzero at a time. The 
point (0,0) corresponds to the SM point and the (*) represents the best fit point. The green region 
corresponds to the 68 percent confidence interval and the red region to the 95 percent confidence 
interval, respectively. The best fit point is doubly degenerate up to a sign flip of the best fit point 
coordinates. 
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(a) 


(b) 


(c) 


Figure 4: Marginalized global fits of the CP-odd parameters with fw = /b = 0. 
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Figure 5: Marginalized global fits of the CP-odd parameters with fw = fsw = 0. 
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Figure 6: Marginalized global fits of the CP-odd parameters with fBB = fsw = 0. 


15 











• We have five CP-odd parameters {fw^ fsB^ Ibw^ fww) ^^nd in the first case we consider 
any two of them to be non-zero and put limits on them. There are a total of ten such 
combinations. The allowed parameter space with 68 % and 95% confidence levels are shown 
in FigJ^ In generating these plots we have varied parameters freely. Even if one considers 
the perturbativity argument, the tightest upper bound on these CP-odd parameters is of the 
order ^ 100. We can see that in all cases the allowed parameter space is bounded. The 
parameters which enter the 77/1 vertex ( 77/1 family) i.e., fBB^fww smd fsw sire in general 
more constrained than those which do not enter the 77/1 vertex (non- 77 /i family) i.e., fs 
and fw- When any of the three 77/1 parameters is taken together with fs or fw as shown in 
Figs [3^3c|3d|3e , [3h|and[3ll we find that it is much tightly constrained and allowed values are 
of 0{1). However, among themselves these parameters are highly correlated and cancellation 
among them leads to larger allowed values of (9(10) (see Figs. [3fpg| and [3j]) . The nature of 
slope in these figures is related to the relative sign among /bb^ fww and /bw- We also note 
that out of fB and fw> fw is always more constrained. For example, the maximum allowed 
value for is ^ 100 while the allowed values of fw are less than 60, see Fig. 3b, This 


observation can be attributed to the relative size of their coefficients in various observables. 
We have also found that the inclination of the plot in Fig. is governed by the C^zh which 
enters Zh and VBF production channels and affects the Stotal [Fq.(| 


• In the second case we consider three parameters at a time in the global analyses. Once again 
there are ten such combinations. Following the general conclusions of two-parameter case, 
we can categorize these combinations into three groups. This categorization is based on the 
number of parameters from the 77/1 family being present in each combination. Thus we 
have one combination where all the parameters are from 77/1 family (Gl), six combinations 
where two parameters are from 77/1 family and one from non- 77 /i family (G 2 ) and three 
combinations where one is from 77/1 family and the other two are and fw (G3). We present 
results for three representative combinations (one from each group): (i) { fBB^fww^fBW 
}; (ii) { JbbJwwJb }; (hi) { JwwJbJw }• The allowed parameter space for these 
combinations are shown in Figs. Hi and respectively. The two-parameter plots here are 

obtained after marginalizing over the third parameter. 


As compared to two-parameter plots (see Fig. |^, the allowed region here is more diffused 
due to the presence of the third parameter. This is particularly noticeable for the first set 
of plots in Fig.J^ We find that fww being present in all VVh coupling gets stronger 
bounds, while fBW is least constrained in set (i). From the plots of Fig. which belong 
to set (ii) we can infer that the parameters of 77/1 family are more constrained. However, 
mutual cancellation still allows values of (9(30) for them. The opposite inclinations of the 
plots in 


5a 


and 


5b 


can be related to how /bb and fww enter in 77/1 vertex to maintain 
cancellation and it is confirmed by the nature of the slope in Fig. Note that Fig. [5^ is 
very much similar to the corresponding plot in Fig. This is expected because the global 
analysis is not very sensitive to /b- Interestingly, the parameter sets {fww/fBw/fw} and 
{fBB^ fBW^ fB} which would also belong to group G 2 are less constrained. We find that the 
parameter regions are still bounded but boundary values for fw and /b in these sets are 
about 800 and 2500, respectively. For the parameters in set (iii) we can conclude that fww 
being the only parameter present from the 77/1 family is very tightly constrained as shown 


by the plots in Figs. 6 b and 6 c 


We have also checked that the constraints on CP-odd parameters from the direct and indirect 


measurements of the Higgs total width [TO, 73 are weaker than those obtained from the global 
analysis. On the whole, / bb , fww and fBW contribute to the 77/1 vertex at tree level and thus 
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are constrained rather tightly compared to fs and fw- The fact that one has SM contribution at 
the one-loop level only is responsible for this, fs and fw are relatively loosely constrained due to 
the lack of sufficient data on the channel h 7 Z, to which they contribute. 

The constraints on the parameter space of our CP-violating operators, as obtained from global 
fits of the (7 -h 8 ) TeV data, are expected to be improved in the high energy runs. A tentative 
estimate of such improvements, as also that in a linear e+e“ collider, can be found, for example. 


in 74 . Going by the estimates of 74 , the uncertainty in the signal strength measurements in the 


next run can be reduced to 33% of the present uncertainties as obtained by both ATLAS and CMS 
for the 77 and WW final states in the gluon fusion channel. A precise answer on the improvement 
of our limits, however, depends also on any possible shift in the central values of the measured 
signal strengths in different channels. This in turn is also a function of the various systematic 
uncertainties in the new run, and therefore we have to wait for more data before some precise 
conclusions can be drawn. A similar consideration applies to a linear collider; precision in coupling 
measurements down to 1 % is expected there in principle (^, but the available statistics as well as 
the systematics need to be known before concrete estimates emerge. 


5 Constraints from EDMs 


The fermionic electric dipole moment receive an additional contribution from these new CP-odd 
higher-dimensional operators involving Higgs and pair of gauge bosons. Nonobservation of any 
fermionic EDMs puts severe constraints on the parameters /^s. The fermion EDM operator is 
defined as. 


- - d/ V’fe) 'ipipi) 


(41) 


where, = §[ 7 ^, 7 ^] and df is known as the fermion EDM form factor. Nonvanishing EDMs 
provide clear hint of CP-violation [75t[7^ . In the standard model, CP violation occurs due to quark 
mixing and it is quite weak (1 part in 1000) [^. On top of that the first nonzero contribution to 
EDM operator in the SM appears at three loop level in quark sector, while, for leptons it arises 
at four loop level. The present upper limits on electron and neutron EDMs are much larger than 
the values predicted by the SM [^. In presence of CP-odd gauge-Higgs operators 77 / 1 , 'jZh and 
WW^ couplings are modified, because of which the leading contribution to fermion EDMs appears 
at one-loop level. Due to this the fermion EDM measurements can provide stringent bounds on 
the CP-odd parameters. Note that the contribution to the fermion EDMs from CP-odd WWh 
and ZZh vertices can result only at two-loop level. Since the two-loop effects are expected to be 
subdominant, we will derive constraints on the CP-odd parameters from one-loop fermion EDM 
calculations. The diagrams that contribute to the fermion EDM at one-loop level are shown in 
FigHI 

The expression for the fermion EDM form factor df at one-loop due to the 77 / 1 , ^Zh and WW^ 
vertices is given by the following equatioi][^ 


where. 


df = 


ai = - 


rrifea 


TTV^ 


diKi{A, rrih) + a2i^2(A, m^, rrih) + mw) 


Qf 

45 ^ 


a 


'j'jhj ^2 


(h/ Qf^w) 


twS 


a 


'jZh') ^3 


2W 


4 


Cww 


7- 


(42) 


(43) 


^We have observed a relative sign change in the contribution of the WW^ diagram to the EDM, for the u and d 
quarks, which is taken care of by the factor If in as. It was unaccounted by Ref. where a similar calculation is 
done. 
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Figure 7: One-loop diagrams contributing to fermion EDMs. The blobs show the effective vertices 
arising out of the CP-odd operators. 


Here, v is the electroweak symmetry breaking scale, a is the fine structure constant, //is the 
third component of the fermion Isospin and Q/ is its electric charge quantum number. We have 
neglected the fermion masses (m/) with respect to other mass scales in the loop. The one-loop 
factors Ki and K 2 are calculated in dimensional regularization (d = 4 — 2e). Since these loops are 
UV divergent, we renormalize them in MS scheme and identify the renormalization scale with the 
cutoff A. The expressions for Ki and K 2 are given by. 


//i(A,x) 


K2{A,x,y) 


V 


1 h? 3 
-In^ + 




4J’ 


1 


x‘^ — 


3i 

(44) 

+ 4 ]- 

(45) 


In the above, the finite factors of | are artifact of dimensional regnlarization. These factors do not 
appear in naive cutoff regularization. 

The latest experimental bounds on the electron and neutron EDMs are (78}|80] , 

\de\ < 8.7 X 10“^® e cm 

\dn\ < 2.9 X 10“2® e cm. (46) 


Note that the EDM contribution is 0{^) in the cutoff scale, therefore, it is expected to provide 
stronger bounds on CP-odd parameters than those obtained from EWP and LHC data. Like the 
electroweak precision observables calculated in section 3, the EDMs also have explicit dependence 
on A. Due to this, we provide EDM constraint equations for three different choices of cutoff scale 
A = 1, 5 and 10 TeV. 


• A = 1 TeV 

|de| = |233.86 /b + 260.45/w - 390.92Jbb - 337.72 fww + 858.63/^^1 < 1 
\dn\ = |7.02/b + 13.8l/m — 8.91/bb + A.OOfww + 9,2.90fBw\ < 1 (47) 

• A = 5 TeV 

\de\ = |13.87/b + 15.63/m - 24.60 /bb - 21.08fww + 52.34/^^! < 1 

\dn\ = |0.40/b + 0.85/m “ 0.57 /bb + 0.33fww + L36/Bm| < 1 (48) 

• A = 10 TeV 

\de\ = |3.95/b + 4 . 47 /m — 7.11 /bb — 6.08/mm + 15.02/Bm| < 1 

\dn\ = |0.11/b + 0.24/m - 0.16 /bb + O.lO/mm + 0 . 39 /Bml < 1 (49) 
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is calculated in terms of constituent quark EDMs using the 

In calculating the above constraints 


81 


The neutron EDM form factor (dn 
relation dn = — ^du^ from the chiral-quark model 

on the EDMs we take a = 1/137 and Mh = 126 GeV. Because of a stronger experimental limit 
on electron EDM, the coefficients of parameters in de are larger than those in Our constraint 
equations for electron EDM form factor differ by an order of magnitude from those obtained in 
Ref. mainly because we have used the most updated experimental bound on electron EDM [78] . 

We first consider the case when any two of the five parameters are kept nonzero. In this case, we 
freely vary the parameters. The constraints on the parameters are shown in Fig for A = 1, 5 and 
10 TeV respectively. In all combinations we get bounded regions. The inclinations of the constraint 
regions can be understood from the relative sign between the parameters in the expressions for the 
EDMs. As expected the constraints for A = 1 TeV are tighter than those for A = 5 and 10 TeV 
and this is corroborated by the size of the coefficients entering in EDM expressions. In A = 1 TeV 
case, the allowed values for parameters can reach (9(1) values at maximum. As we push the cutoff 
scale higher, the allowed range for parameters also increases. For example, for A = 10 TeV the 
allowed values can become (9(10) or larger in some cases. 

A naive comparison with the EDM calculation carried out in Ref. with only 77/1 CP-odd 
coupling suggests that for A = 1 TeV, the constraint equation for electron EDM measurement 
would be \fBB-\- fww\ ^ 0.0036. In the presence of CP-odd 7 Z /1 coupling, which arises from same 
operators, this constraint equation would change. For example, our calculation for electron EDM 
constraint implies, I/bb + 0-86fww\ ^ 0.0026. In both the cases, the parameters would be allowed 
to take very large but fine-tuned values. However, we would like to point out that after including 
the constraints from the neutron EDM measurement, these parameters cannot take values larger 
than (9(0.1). 

When we take three parameters nonzero at a time, the parameters are scanned in the range -200 
to 200. In Fig[^ we give two dimensional projection plots of the three dimensional constraint region 
for A = 1 TeV. For comparison purpose, we present the plots in three categories discussed in the 
global analysis. We can see that with three parameters present, the constraints are more relaxed 
than when only two of them are nonzero. However, the parameters are still quite correlated, /b 
and /bb can often reach the boundary of the scanned regions. In fact it is very difficult to obtain 
closed boundaries. As we increase the range for parameter scan the allowed values for CP-odd 
parameters become very large [(9(1000)]. However, it is important to note that for too large values 
of parameters, the two-loop EDM constraints may become relevant and, therefore, should also be 
taken into account. As we turn on more parameters, the correlation among parameters constrained 
by EDM is relaxed and the allowed parameter space also expands. 

From the experimental perspective, a number of new EDM experiments promise to improve the 
level of sensitivity by one to two orders of magnitude in the coming years. For example. The 
Institut Laue Langevin (ILL) cryogenic experiment and the Spallation Neutron Source (SNS) 
nEDM experiment ( [ 8 ^, 


, [84] ) aim at improving the upper limit on neutron EDM by two 
orders of magnitude, i.e down to (9(10“^^)e.cm. This would imply that the numerical coefficients 
in the constraint equations for dn in Eqs. ( [48| ) and ( [4^ would become stronger by almost 

two orders of magnitude and thus the allowed parameter space for f's will be even more severely 
constrained, unless, of course, there is direct evidence of neutron EDM in the aforesaid experiments. 


6 Discussion 

We now highlight the important features of our analysis presented in sections 3, 4 and 5. We try 
to draw a comparative picture and address some of the issues relevant to the analysis. 
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Figure 8: EDM constraints keeping two parameters nonzero at a time for three representative 
values of A = 1 (blue), 5(green) & lO(red) TeV. 
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(a) 





(b) 





(c) 

Figure 9: EDM constraints with {&)fw = /s = 0, {h)fw = fsw = 0 and (c)/ss = fsw = 0 for 
A = 1 TeV. Parameters are varied in the range between -200 to 200. 
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• In the two parameter case, among all the constraints, the EWP constraints on CP-odd 
parameters are always the weakest while EDM constraints are the strongest. In general, the 
correlation among parameters is stronger in EWP and EDM cases as compared to the LHC 
case. 

• In case of global fit of LHC data whatever contributes to /i ^ 77 receives stronger constraints, 
since here the “tree level” contributions from the dimension six operators are essentially at 
the same level as the “one-loop” SM contributions. There are three such parameters namely 
Ibb^ fww 2 md /bw- Only a strong cancellation or correlation can keep the individual values 
of these parameters large. 

• The same should apply to parameters contributing io h ^ 7 Z channel which also includes 
fs and fw- However the limits based on current data are rather weak there and we have 
not included them in our analysis. The global analysis puts some limits on fs and fw as 
the yZ/i coupling modifies the Zh and VBF production channels and the total Higgs decay 
width. 

• The channels dependent on the WWh and ZZh couplings yield relatively weak constraints 
from the global fits as a whole, since the higher dimensional interactions are inadequate to 
override the tree level SM contributions. 

• In two parameter case, LHC data bounds on /bb, /bw and fww are stronger when any of 
these is combined with /b or fw- These bounds are comparable to corresponding bounds 
obtained from EDM measurements. 

• In three parameter case, the bounds from LHC data are stronger than those obtained from 
EDMs. In fact, for some parameters (for example, /b) values of the order 1000 are also 
allowed by EDM data. However, the parameters constrained from EDMs still display a tight 
correlation. 

• If we allow dimension -6 CP-even operators to coexist with the CP-odd ones, the EWP and 
LHC observables would receive contributions at 0(1/A^) (as a result of interference with 
the SM) as well as at 0(1/A^). Since there is no interference between CP-odd and CP-even 
operators in total rates, the CP-odd interactions studied by us should contribute to signal 
strengths on the order of ( 1 /A^) . For the sake of consistency, at (9(1/A^) the contribution 
from dimension -8 CP-even operators via its interference with the SM should also be consid¬ 
ered. Thus, in the presence of CP-even operators, the constraints from LHC and precision 
data can be relaxed. However, under the assumption that the effect of CP-conserving new 
physics is not significantly large, the constraints obtained by us on CP-odd operators are in 
a way the most conservative estimates of the allowed parameter space. Of course, the limits 
from EDM analyses remain the same in all the cases. 

• A comparison between the relative strengths of the EDM and LHC constraints, we emphasize, 
is most transparent only when the CP-violating operators alone are considered, since the CP- 
conserving ones have no role in EDMs. Their inclusion, albeit via marginalization in the LHC 
global fits, will serve to relax the corresponding limits beyond what we have obtained. But 
then, it ceases to be a one-to-one comparison between the two kinds of constraints, something 
that we have intended to do from the beginning. 

We have already seen that the Lorentz structure of anomalous CP-odd couplings is unique 
and these are just some linear combinations of the CP-odd parameters (see Table [^. Since in all 
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Couplings 

LHC data 

EDM 

2P case 

3P case 

2P case 

3P case 

\CwWh\ 

0-60 

0-60 

0 - 0.17 

0-55 

\Czzh\ 

25 - 100 

25 - 80 

0.11 - 0.20 

0.15 - 33 

1 I 

0 - 0.8 

0 - 0.5 

0 - 0.16 

0.02 - 52 

\C^zh\ 

20 - 25 

15 - 25 

0.03 - 0.25 

0.05 - no 

\CwW'i\ 

0-40 

15 - 40 

0 - 0.15 

0.02 - 47 


Table 4: Limits on CP-odd coupling strengths from LHC data and EDM measurements for A=1 
TeV. 2P and 3P stand for two parameter nonzero and three parameter nonzero cases respectively. 


the observables these couplings enter directly, it is instructive to know the kind of values these 
couplings can take as a result of our analysis presented above. In the calculation of 5, T and U 
parameters all the CP-odd couplings directly enter. In the global analysis only Cyvh couplings 
participate, therefore, constraints on Cwwv from the LHC data are indirect. In EDM calculations 
only Cjjh^CjZh si^rid Cww^y enter directly and therefore limits obtained on Cwwh smd Czzh sire 
also indirect. The limits on the strengths of the anomalous couplings are listed in Table Since 
electroweak precision bounds on fi are the weakest in the table we compare bounds on couplings 
due to LHC data and EDMs. The comparison is presented for both the two parameter (2P) and 
three parameter (3P) nonzero cases. Since Cwwz is proportional to CwW'y^ fhe limits on Cwwj 
can be easily translated to limits on Cwwz- Looking at the LHC limits on the couplings we 
find that is the most constrained coupling. Also, the LHC limits in 2P and 3P cases are 

comparable. On the other hand the EDM limits on couplings in 2P case is always stronger than 
in 3P case. In 3P case the lower limits result from the three parameter sets {fB^fsB^/Bw} ^^nd 
{fw^ fww^ fBw}- The parameters of these sets are found to be very fine tuned. It is important 
to recall that in 3P EDM case, the parameters are varied in the range between -200 to 200. We 
find that as we increase this range, the maximum allowed values for couplings also increase. We 
would also like to point out that the correlations among CP-odd couplings are mostly similar to 
those found among CP-odd parameters. 

The triple gauge boson couplings (TGCs) WW'j and WWZ can also be constrained using 
the collider data on gauge boson pair production. Data from Tevatron and LHC are used mainly 
to constrain the CP-even anomalous couplings as the observables used are not sensitive to CP- 
odd couplings [85}|^ . On the other hand, the experimental analyses at LEP which studied the 
angular distribution of final state particles are sensitive to CP-odd TGCs. A comparison between 

^Although the EWP constraints are weaker in general, we find that the limits on Cwwv from EWP and LHC 
data are comparable in 2P case. 
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the CP-odd sector of TGC Lagrangian 54,91,92 


and our effective Lagrangian (Eq. implies 


CwW'i 

~1^ 


Cwwz 

A2 


Sw 

o 

m 




w 

cw ~ 


m 


w 


(50) 

(51) 


Using the relation between Cww-f stnd Cwwz we get, kz — 


LEP limits on kz are 79 


k^. At 68 % CL, the combined 


-0.14 <kz < -0.06. 
These limits when translated on Cwwz and CwW'y become. 


(52) 


-0.19 < 
0.15 < 


[TeV^] < -0.08, 
[TeV“2] < 0.36. 


(53) 


Note that these limits are comparable to the limits obtained from EDMs in 2P case, however, 
information on the sign of the couplings is also available. 

Other than the VVh vertices considered in our analysis, quartic VVhh vertices also arise out of 
gauge invariant CP violating operators. One can thus expect some correlated phenomenology from 
the trilinear and quartic interactions, since the former arise essentially on replacing one Higgs by 
its vacuum expectation value in the quartic terms. In the analysis presented in this work, the 
production and decay channels considered are not affected at tree level by such quartic VVhh 
vertices. Thus in the context of the present analysis, any constraints on such quartic vertices from 
observed data are likely to be weaker than those obtained from the VVh{V — effective 

interactions, since the quartic couplings would entail Higgs pair production. For example, the 
VVhh vertex may contribute in addition to the hhh vertex toward a di-Higgs final state. However, 
one needs to wait for a large volume of data on Higgs pair production to see such correlated 
phenomena. In general, limits stronger than what we have obtained are not expected. Also, the 
contributions from such VVhh CP-odd vertices to EDMs and EWP observables come at higher 
loop levels and thus are expected to be substantially weaker than what we have obtained for the 
trilinear terms. 

It is a well-known fact that the observed baryon asymmetry in our universe cannot be explained 
by just the CP-violating phase of the Cabibbo-Kobayashi-Maskawa (CKM) matrix in the SM. 
The presence of additional sources of CP-violating operators arising from the anomalous VVh 
interactions may in principle explain the observed baryon asymmetry of the universe. However, 
a more careful scrutiny of this picture reveals that our CP-violating operators are not sufficient 
to trigger strongly the first order electroweak phase transition required for the baryogenesis. For 
this, one has to extend the Higgs sector of the SM by introducing new particles which couple to 
the Higgs boson and thus modify the Higgs potential such that it leads to a strongly first order 
electroweak phase transition [^. 


7 Summary and Conclusions 

We have analyzed CP-odd VVhiV = VE, Z, 7 ) and WWViV = Z, 7 ) interactions in terms of 
gauge invariant dimension -6 operators, obtained as the artifacts of physics beyond the standard 
model. The most complete set, comprising five gauge-Higgs operators, has been taken into account. 
We have derived constraints on the coefficients of such operators using electroweak precision data. 
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LHC data on Higgs and limits on the electric dipole moments of the neutron and the electron. 
With A as the scale suppressing the CP-violating operators, precision parameters as well as LHC 
observables receive contributions ^ 1 /A^ from the CP-odd couplings, while contributions ~ 1 /A^ 
to EDM observables are expected. The constraints obtained from the T and U parameters are 
the weakest, while the bounds from EDMs are the strongest with two nonvanishing operators. The 
global analysis of Higgs data from the LHC puts stronger constraints on those CP-odd effective 
couplings which contribute to /i ^ 77 , as compared to those which do not. We also indicate 
situations where large values of certain couplings are allowed by all constraints, when they appear 
in combination. The constraints coming from LEP on CP-odd form factors CwW'y and Cwwz are 
consistent with our limits obtained from EDMs in the case when any two out of five parameters are 
nonzero. It may be of interest to find out new physics scenarios where, by integrating out heavy 
degrees of freedom, one may arrive at large correlated values of such operators. In a subsequent 
work, we hope to discuss some observables that may help one in probing these operators in 
the 13 and 14 TeV LHC runs. 
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Appendix A Gauge boson two-point functions in presence of CP- 
odd couplings 

We define the two-point function for electroweak gauge bosons Vi and V 2 as, 

= g'^^Uv.vAP^) +p^p^Uv,vAp^). (54) 

In electroweak precision observables only IIy^v 2 {p^) contribute. Below we give their expressions 
due to CP-odd couplings in terms of one-loop scalar functions Aq and Bq. 




= (Sy {C^jh -y) (2m^ +y) 

+ C^zhMmz) + - ^C^jh'mhBo{p^,0,ml) - SC^zh (m| +/) 

+ m Bo{p^,ml,ml) + SAo{ml) (7^^ (ml+p^) + C^zh W “ "^1 +P^)) 

+ ^C^zh^oiml) (m| - ml) + p^ {7Cl^h {P^ “ + TC^zh {P^ “ 3 + m|)) 

+ 2Cl^^^ (l2m^ + 7^))) , (55) 
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2 2 

(-^C^'fhCjZh {ml - p^f Bq{p^, 0, ml) - SCyZhCzzh {ml - 2ml (m| + 

+ {ml -P^Y) Bo{p^,m%,ml) - QCwwyCwwzp'^ (2m^ +p^) 

+ ^C^ZhAoiml) {C^^h {ml+p^) + Czzh {ml - ml+p^)) + SCjZhCzzhAoiml) 

{—ml + ml +p^) + l2CwW'rCwwzp‘^Ao{ml^) + p"^ {7p‘^{CjZh{Cjjh + Czzh) 

+ ‘2‘Cww-iCwwz) — ‘2‘IC^zh {mliC^^h + Czzh) + Czzhml) + 24CwwjCwwzmly)) , 

(56) 


B-zziP^^) — (3p^ {C^zh Bo{p‘^,0,ml) — 2Cl/wz i^mlr +p^) BQ{p^^ml^,mlr) 

+ + Clzh^o{ml)) — 3C^zh''^h^oip‘^,0,ml) — SClzh {^1 ~ ‘^ml {ml +p^) 

+ {ml-p^f') Bo{p^,ml,ml) + SAo{ml) {C^zh {^l+P^) + C'lzft {ml - ml+p^)) 

+ ^Czzh^o{ml) {ml — ml) +p‘^ il^Clzh {p^ ~ ^ml) + 2C^r^rz (l2m^ + 7p^) 

+ '^Clzh {P^ ~ 3 {ml + ml)))) , (57) 


2 2 

niyiy(p^) = (^~^wW'r “7*^) {mw +7*^) -Bo(7'^,m'^,0) + C'^|4/;j(—p^) {ml 

~ ‘^ml {mil + p^) + {m^Y — p^) ^ Bo{p‘^, m^y, ml) — Cliyiz {mw ~ mw {^ml + p^) 

+ mw {ml + Smlp^ - (p^)^) + p^ {ml - Bo(p^m^,m|) - Cwwz^oiml) {mw + p^) 

{^w ~ ~ + 3Ao(m^) {j^w ~ +P^) 

+ C^yyz {'^W ~ '^W {'^‘z + lOp^) {p^ — '^z))) + ^^WWhP^^oi'fT^h) {'^h ~ '^w P^) 

+ {Cww-f ~ + 7(^^)^) + IC^yy^p^ {p^ — 3 (m^ + rn^)) + C^yyz 

(-7p^ (2m^ + 3m|) + 87m^ (m^ + m|) + 7((p^)^))) . (58) 


Out of these, and lizz vanish at p^ = 0. Note that in liww there is an overall l/p‘^ 

dependence. We would like to mention that both I\^ww and its derivative converge smoothly in 
^ 0 limit. The one-loop scalar functions in n = 4 — 2e dimensions are given by, 


where. 


AQ{ml) 


Bo{p^,ml,m\) 



IGtt^ 


ml 


I 


(Bl 

(27r)" 


1 ri 

IGtt^ e 


1 

P — mg 

t + 1 - In(mo) , 

1 

{B -ml) {{I +p)2 
- A{p‘^,ml,ml) , 


mt 


A{p^,ml,ml) 



In —x{l — x)p^ + x{m\ 


ml) + m^] . 


(59) 


(60) 


(61) 
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This form is suitable for computing Sq and its derivative with respect to at = 0, which we 
require to calculate 5, T and U parameters discussed in section 3. 
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